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Abstract 

Prediction of aerodynamic loads on bodies in 
arbitrary motion is considered from an acoustic 
point of view, i.e., in a frame of reference fixed 
in the undisturbed medium. An inhomogeneous wave 
equation which governs the disturbance pressure is 
constructed and solved formally using generalized 
function theory. When the observer is located on 
the moving body surface there results a singular 
linear integral equation for surface pressure. Two 
different methods for obtaining such equations are 
discussed. Both steady and unsteady aerodynamic 
calculations are considered. Two examples are 
presented, the more important being an application 
to propeller aerodynamics* Of particular interest 
for numerical applications is the analytical 
behavior of the kernel functions in the various 
integral equations. 

Introduction 

Aerodynamics is a fascinatinq but somewhat 
strange field. It is fascinating because it con- 
tinues to offer many interesting practical problems 
which must be solved or are in need of better solu- 
tions. It is strange compared to other fields of 
science because most available analytical methods 
of attacking aerodynamic problems are indirect in 
the following sense. In aerodynamics one rarely 
starts with a differential equation and boundary 
conditions and then proceeds to solve a problem. 
Of course this state of affairs occurs because no 
general analytical method of solution of the Euler 
or Navier-Stokes equations is known* Even when 
linearized theories are used, quite often the 
direct method of solution is avoided. Instead, one 
starts with doublet or vortex distributions over a 
surface and derives an integral equation in terms 
of the unknown source strength. These methods, 
although correct, are somewhat unsatisfactory be- 
cause they require a great deal of physical insight 
and knowledge of physics. As a result, engineers 
and students often find aerodynamics a much more 
difficult subject than need he. It is one of the 
purposes of this paper to propose a method to rem- 
edy this circumstance through the use of acoustic 
equations. We limit ourselves to linear aerody- 
namics only. Our main objective is to obtain re- 
sults for prediction of aerodynamic loads on rota- 
ting blades. A comprehensive review article on 
this subject has been written recently by Johnson 
[13* 


It has long been known that linear aerody- 
namics and acoustics are closely related [23. This 
fact has been exploited before by many researchers 
such as Garrick [3] and Kondo [41. More recently, 
Hanson has specialized his acoustic formulations 
for rotating blades to apply to aerodynamic predic- 
tion [5] . His velocity potential approach is 
closely related to that of Kondo [4J and involves 
infinite series of special functions which can 
cause computational problems. It is also restric- 
ted to propeller-like motion and thus cannot be 
adapted easily to propellers with non-uniform in- 
flow or to helicopter rotors. 

Das has made significant contributions to the 
application of concepts from acoustics to linear 
aerodynamics [6]. His work is closely related to 
the first of two methods presented here. His 
approach, however, is from the point of view of 
singularity distributions which for reasons discus- 
sed above is avoided here. The approach presented 
here also bypasses much of the lengthy algebraic 
manipulation common in singularity distribution 
methods. 

The approach we propose here is based on the 
solution of the wave equation in the time domain 
and in a frame fixed to the undisturbed medium. 
Consequently, compressibility effects are automat- 
ically included in the analysis. This unconven- 
tional approach can cause some difficulty for those 
who are accustomed to working in a co-moving frame, 
i.e., a frame fixed to the moving body. However, 
this difficulty is compensated for by results of 
great generality and usefulness. Starting with an 
acoustic equation for a body in arbitrary motion, 
we can derive an aerodynamic integral equation by 
moving the observer onto the surface of the body. 
Both steady and unsteady motion are considered. 
Two different methods for deriving an integral 
equation are discussed. Surprisingly, the simpler 
one applies only to unsteady periodic and steady 
camber problems but not to the steady angle of 
attack problem. Essentially, the unsteady aerody- 
namic problem is an acoustic problem in a sense 
that will be explained later. 

The origin of this paper is a set of lecture 
notes written for a short course on acoustics and 
aerodynamics of propellers at the von Karman Insti- 
tute [73 . Long has demonstrated the usefulness of 
the acoustic approach by deriving an integral equa- 
tion and working out several numerical examples 
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[8] . More recently Milliken has followed a related 
approach by deriving the aerodynamic kernel func- 
tion for a rectangular panel in steady rectilinear 
motion using the acoustic method [9). This result 
will be presented here-- as preparation for a more 
complicated kernel function for propellers. 

Following a discussion of the governing equa- 
tion of acoustics and its relevence to aerody- 
namics, application to aerodynamic problems will be 
considered* Two examples will be given, the second 
beinq an application to propeller aerodynamics. 
For the case of the propeller, an analytical treat- 
ment of the singularities appearing in the integral 
equation is given. As a result, only non -singular 
integrals are required for numerical work. 


the time derivative of both sides of Bq. (5) and 
using Eqs . (2) to (4), we obtain 

□ 2 p=V* [Apn|vf|«(f )] (6) 

Here n=7f/|vf| is the unit normal to f=0 pointing 
toward the suction side. Equation (6) is a special 
case of the Ffowcs Williams-Hawkings (FW-H) equa- 
tion [7, 10, 13]. It is the governing equation for 
all acceleration potential (doublet lattice) meth- 
ods. To relate the above results to those of clas- 
sical aerodynamics, consult reference 7, The math- 
ematics of this section is discussed in references 
10 to 12. 

The Acoustic Approach To Aerodynamics 


The Governing Equation 

Let <f> (x, t) and p (x, t) be the velocity .poten- 
tial and the pressure, respectively. Here the x- 
frame is fixed to the undisturbed fluid medium in 
which the speed of sound is c. The velocity poten- 
tial in the linearized form satisfies the wave 
equation 

!_ JL| . v 2 * z nVo (D 

c at + 

The fluid perturbation velocity is u=V$ , and the 
relation between p an $ is described by 


We consider two different methods in this 
section. The first leads to a rather conventional 
integral equation. We begin with the formal solu- 
tion of Eq. (6), which is, 

4np(x,t)-9* / -1 Apn|vf |i(f)«(g)dydT ' (7) 

where g»r -t+r/c, r=| r | = | x-y | , and (x,t) and 

(y, t) are the observer and source space-time vari- 
ables, respectively. For aerodynamic application 
we use one special form of the various equivalent 
interpretations of Eq. (7) [14). Taking the diver- 
gence term in Eq. (7) inside the integral and using 
the relation 


P=-P 0 


3t 


where p is the density of undisturbed medium, 
o 


(2) 


3 . 6(g) . ! 3 A*** 

3x t L r J " c 3t l r r 2 


( 8 ) 


Consider a lifting thin body in arbitrary 
motion. In the following we will omit consider- 
ation of the thickness effect. Within linearized 
theory, boundary conditions on the body are re- 
ferred to an infinitesimally thin mean chord sur- 
face whose velocity is everywhere tangent to it- 
self. Assume that f(x ,t)=0 describes this surface 
in motion such that Vf points toward the suction 
side of the lifting body. Let k(x ,t)=0 describe 
the wake of the body, also assumed infinitesimally 
thin. Across the surface f=0, the jump in surface 
pressure, Ap=p lower ~P upper' is related to the 
velocity potential by Eq. (2) as follows 

Ap-p 0 A* t =0 (3) 


where r^r^/r ,we get 

4 irp ( x , t ) = - 1 | 5(f )5( g )dydT 

C 3t r ( 9 ) 

- /*P|£li|v£| 6 ( £ ) S ( g )dydT 
r 

The interpretation of the integrals in Eq. (9) 
is given elsewhere [14], Let 2 be the surface 
generated by the intersection of f«0 and the col- 
lapsing sphere g=T-t+r/c=0. Since f=0 is assumed 
to move tangent to itself, M n =0 on the 2- surface 
(therefore A in reference 14 is equal to 1) and we 
can write Eg. (9) as follows: 


where $ t -3<fr/3t. Across tlv~ wake, since Ap«0, erne 
gets 


4irp(x, t )= - ~ |^/ i[Apcos8) T ^d2 


F=0 


( 10 ) 


A<J> t =0 (4) 

Note that ^Vupper’* Viewer in Eqs. (3) and 

Following reasoning used in classical aerody- 
namics, has a jump across f=0 and k=0. We wish 
to obtain a wave equation in which this discontinu- 
ity is included explicitly in an inhomogeneous 
source term. To do this we replace the derivatives 
in Eq. (1) by generalized derivatives [7, 10-12). 

The result is 

□ 2 <fr=-V* [A* Vf <S( f ) ] -V- [ A* Vk«(k)] (5) 

where the bar over the operations denotes general- 
ized differentiation and 6(* ) is the Dirac delta 
function. Here again - <f> lower * Taking 


- / ^j(Apcos0) T# d2 
F=0 r 

where F (y; x, t ) = [ f (y, t) ] ret =f ( t-r/c) and T*is the 
emission time. One more step is needed to get the 
required result. Using Eq. (2) on the left side of 
Eq. (10) and then integrating both sides with re- 
spect to the observer time from - • to t, we obtain 

4irp A(x,t)= — / — (ApcosQ) d2 

° ° F=0 r T (11) 

+ / / -^-tApcosej^dldt' 

— F=0 r 

where T** is the emission time corresponding 
to (x,t'),t' being the integration variable with 
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respect to observer time. Upon calculating the 
normal derivative 3/3n of both sides of Eq. (11), 
we qet the aerodynamic integral equation 

4l ' P 0 V n (5 ' ,t) “ c In-f 7(AP cos0 I T . dj: 

0 n c 3n F =0 r T (12) 

a t , 

+ -37 // -jtApcosej^dSdt' 

F=0 r 

where v is the local normal velocity of the sur- 
face. ^e right side of Eg. (12) must be inter- 
preted as a limit for x approaching the surface 
f=0. 

Equation (12) is a linear singular integral 
equation. It is the basis of the acceleration 

potential or doublet lattice method. We note, 
however, that there is no restriction made in Eq, 
(12) to uniform rectilinear motion of the surface 
f=0. The aerodynamic kernel is derived from this 
result and will be discussed later. We mention 
here that the second integral on the right of Eq. 
(12) is an improper integral. Thus the normal 
derivative cannot be taken inside without careful 
analysis since it leads to a singular integral. 


kernel function. This can be achieved most effi- 
ciently using some elementary generalized function 
theory. 

^ Let us define the following integral in which 
q (y;x,t) is an arbitrary function and F^ has the 
same meaning as defined above: 

i 

Note the surface integration is with respect to the 
variable y and the normal derivative is with re- 
spect to x , On the surface f=0, let the ith panel 
be specified by f^-0, £<y,t)>0 where f^=0, l =0 
specifies the edge of the panel. Then using L 

(y;x . , t )=£(y, t- — ) ,r=|x . -yi , Eq. (14) can be 

3 c * J 1 

written as 

i= -§7f H(L)c ' d5: 0 5) 

where H(L) is the Heaviside function. The limits 
of the integrals in Eq. (15) can be considered 
constants now and the normal derivative can be 
taken inside as follows: 


A second method of using acoustic equations in 
aerodynamics is to treat Eg. (10) itself as an 
integral equation on the unknown surface pressure 
as suggested by Long [15] and Farassat [16]. For 
reasons which will be explained in the next sec- 
tion, this method fails for the steady angle of 
attack problem. One is thus obliged to use the 
first method. However, in the case of periodic 
motion the second method appears to lead to a 
useful integral equation whose solution is 
considerably simpler than that of Eq. (12). 

Further Analysis and Applicability to 
Steady and Unsteady Motion 


1 = 5 ( L ) d E+ /H ( L )|2d Z 


1 3n 


|2d 

3n 


(16) 


We also have, from the definition of L, 

_=-_ cos 0_ 07) 

in which 0“" is the angle between n J.at the j-th 

panel) and the unit radiation vector r (from the i- 
th panel to the j-th panel). 


The first integral in Eq. (16) is a line inte- 
gral over the edge of the £- surface generated by 
the i-th panel. It can be written in a compact 
form by first writing Eq. (17) as follows: 


Equation (12) can be used for numerical calcu- 
lations by dividing the wing or blade surface into 
panels and approximating the integrals assuming 
that Ap is uniform on each panel: 


4lrp oV V t)= £ 4 P i { ^ ‘S' F / =0 > os9, T * d5: 

i 

+ I 1 / “4? tcosQl^dldt'} 


(13) 


F. =0 

l 

lisp, K. . 

K 0 l^i l] 


Here F^=F^ (y;x.,t)s*0 is the surface generated by 
the intersection of the i-th panel with the col- 
lapsing sphere g= r -t+r/c-0. The function ■ is 
the aerodynamic kernel function. Physically, it is 
the normal velocity induced at the j-th panel by 
the i-th panel. 


■E =M v Mcos9' 08) 

where M^M* v is the local Mach number in the direc- 
tion of the inward geodesic normal v to the edge of 
the panel f.=fs=o. The geodesic normal v is defined 
as the unit vector which is perpendicular to the 
edge and tangent to the panel f . — Jt=0 • Then, let 
1 2 

(u , u*) be an arbitrary curvilinear coordinate 
system on £ and q^j be the determinant of the 

coefficients of the first fundamental form. We 

r *1 o 

have dj> / du • Using the coordinate trans- 
formation u^+L , we have 


| VL. | A (2) du 1 du 2 = 

M=Ma 0 


IHS 

| 3L/3U 1 ) 


du dL =dYdL 


(19-a) 


( 19-b) 


Under some conditions, the integrals in the 
kernel function can be evaluated analytically. We 
present an example in the next section. Quite 
often, however, one has to resort to numerical 
integration but treat the singularities in the 
integrals analytically. In preparation for this 
step, we carry the analysis of the previous section 
further. First we note that F^ is a function of n 
so that essentially Leibniz rule for differention 
of an integral must be used on the integrals in the 


A 2 0 =M 2 cos 2 \Jh-( 1-Mp* r^sin^) 2 ( 1 9-c) 

Here i? is the projection of the local Mach number 
vector^ on the plane normal to the edge curve, 

Mp= |M p j and rp is the unit vector along the pro- 
jection of r on the plane. Also ij; is the angle 

+ > 

between the radiation vector r and the edge turve 
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of the £- surface. Figure 1 illustrates these 
geometric quantities. In Bq. (19-a), dy is the 
element of length of the edge curve of the £- sur- 
face. Substituting Eg.- (18) into the first term of 
Eg. (16), using Bq. (19) and finally integrating 
with respect to L gives 


Note that the index j is implicit on the right side 
of Eq. (23) by taking x=x . in every term of the 

integrands. An example of the use of the above 
equation to find the kernel function for propeller- 
like motion is given in the next section. 


1= / qM v d y+ /H(L)|2dE (20) 


where we define F?=[ F. (y; x, t ) ] „ . 

i i n=u 

In general the second term in Bq. (20) is 
singular for aerodynamic problems. This would 

happen for Bq . (14) if we assume that the integral 
in Eg. (14) is an improper integral. Mathemati- 
cally, 3/3n could not be taken under the integral 
in the classical sense. However, we can interpret 
3/3n as a generalized derivative. Taking the 

finite part of a divergent integral in many prob- 
lems of aerodynamics is equivalent to 
assuming 3/3n is a generalized derivative. We can 
avoid this route and the associated algebra and 
computational difficulties by taking 3/3n out of 
the integral as follows: 


( 21 ) 


We thus write Eq. (20) in the form 

I= / qM v £2 r : dY + It / qdl (22) 

o v A o 3n 0 „ 

F. =0 F. =0 

L=0 L^0 


Note that L=[L] n _ Q also but we do not have to use a 

new symbol. In practice, the second integral is 
evaluated analytically in the region near where 
singularities exist and the normal derivative is 
then taken analytically. In the other regions, we 
simply evaluate J3q/ 3n d£ numerically. 


Going back to the kernel function defined 
in Eq. (13), we note that both integrals in K^j are 

of the type in Bq. (14). However, the first inte- 
gral is proper (or convergent) while the second is 
improper. Keeping the comments of the above para- 
graphs in mind, we use the analysis of Eq. (14) to 
write the kernel function as 


M cos9'cos0 

4 * p 0 K ij”.c f r [ ^T l* d T 

F. =0 
L 1 =0 


+ 7 l { 7 [cose, T . }d5: 

f i-o 

L>0 

. , M cos 9 'cos 0 

+ / / ~2 A V.-**' 

- F?=0 r ° 

L X =0 


(23) 


-§77 J 1 / tcos0] dEdt' 

„0 _ r 
F. =0 

L^0 


Since our procedure is equivalent to what is 
usually done in acceleration potential methods (16, 
17] it is obvious that it is applicable to both 
steady and unsteady motion. We mention here that 
the time integrations in Eg. (23) must, in general, 
be carried out numerically, which would require 
excessive computation time. It would be highly 
desirable if somehow this could be avoided. We 
study the possibility of removing the time integra- 
tion next. 


As mentioned in the previous section, one may 
treat Eq. (10) as an integral equation on p. Since 
we are interested only in the lifting problem, we 
can assume P upper = -P loW er' A P= 2 Plower* 

Now, it would appear that use of Eq. (10) would be 
an attractive method because it seems to say that, 
except for information at the emission time t*, the 
past history of the motion of the body is irrele- 
vant to the determination of the surface pressure. 
However, Ap(T*) in turn depends on Ap at earlier 
times, and in fact the entire history is required 
for one to apply Bg. (10). For arbitrary motions 
of the body, one infers that only the first method 
leads to a proper integral equation unless it hap- 
pens that Ap has some sort of symmetry property in 
time. The most common situations in which this 
kind of behavior exists are for steady and periodic 
surface pressure. Below we will show that the 
second method also breaks down for the case of 
steady surface pressure due to angle of attack but 
not for the camber and the unsteady (periodic) 
cases • 


Long (8] and Farassat [15] have taken two 
different but equivalent expressions of Eg. (10) in 
which 3/3t is taken inside the first integral ana- 
lytically, and then have treated the result as an 
integral equation on p. These integral equations 
are quite complicated. Here we derive by a direct 
method from Eq. (10) the thin body approximation of 
the results of Long and Farassat. First, however, 
we apply Eq, (22) to the two integrals of Eq. (10) 
which results in 


4 Vn (X+ ' t) 4 [ 7 l 


1 ApM^cos 9 'cos 9 




°_ “0 
£=o° 


+ 7 / { 7 [Apcos9, T* }dj: 

£>V 


(24) 


t , 4pM cos 8 'cos 0 

+ / / V— 4 

-« 0 r A 0 T 

£ =8 


+ / / ^[ApcosO] _d£dt' 

-« 0 r 

£> 0 ° 


£ 
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where- now F°=[F(y; x, t) 3 and L=[ l(y, t- — ) ] , 

n=u c n=u 

2(y, T )=0 describing the blade planform. Again we 
have kept 3/3n outside the last integral since 
this integral is improper. 

We now take 3/3t of both sides of Eg. (24), 
We can write 3v /3t in the frame attached to the 
body as follows: 


modified to obtain the aerodynamic kernel. It 
appears, therefore, that, from a computational 
viewpoint, the steady angle of attack problem for a 
thin wing is more difficult than the camber or the 
unsteady loading problem. For the angle of attack 
problem, one is obliged to use Eq* (12); no alter- 
native exists. 

Examples of Application 


3v 8v n 

n y n 

3t “ n 3a 


(25) 


where v is the rate of change of v for an observ- 
n n 

er fixed on the body, V is the local forward speed 

(of relative wind) and 3/3 a is the directional 

derivative in the direction of forward speed. 

Equation (20) then gives 


We present two examples here for the steady 
state case. The first is the determination of the 
aerodynamic kernel function of a lifting surface in 
uniform rectilinear motion. This function can be 
found analytically for a rectangular panel. The 
second example is for a propeller in uniform axial 
motion. For this case we treat only the more dif- 
ficult question of singularities of the kernel 
function. 


. .. 1 3 t 1 ApM^cos 9'cos 8 

4 * p 0 ( V V W )m cTif r ] x* dy 

E=o° 

+ ^7 [apcos9, T. }dz 

£>V <26) 

„ ApM cosQ'cosQ 

E=o° 


+ 2 2t ApcoseI T* d!: 

„0 r 


£>o° 


Considering the steady state case 

first (v =0) , we note that 3v /3o=0 for the anqle 
n n 

of attack problem so that the left side of Eg. (26) 
is zero. We are led to the conclusion that the Ap 
distribution for this case is an eigenfunction of 
the integral operator on the right side of Eq. (26) 
corresponding to the eigenvalue zero. This fact is 
of little help in finding the pressure distribu- 
tion. However, the integral equation for the cam- 
ber problem is well defined and appears useful 
since V3v /3a * 0. This fact must be explored 
further. 

For the case of the unsteady periodic loading 
problem, we can assume both v^ and Ap are propor- 
tional to exp (ioit), where gj is the frequency of 

oscillation, and once again we have an integral 

equation which is nondegenerate since v tO . This 

n 

integral equation is in terms of the unknown com- 
plex amplitude of Ap . To the authors' knowledge, 

this equation has not been used before. 

The chief advantage of using Eg. (26) instead 
of Eq. (12) is that no time integral appears in the 
former. As a result, Eq. (26) requires much less 
computer time than Eq. (12) in the case of the 
camber and unsteady periodic loading problems. We 
can define a new aerodynamic kernel function based 
on Eg. (26). This is done by using a unit ampli- 
tude oscillation of Ap and finding the induced 

v n at a given point. The calculation of the inte- 

arals in Eg. (26) is very similar to conventional 
acoustic calculations for rotating blades [14] 
which are at present in the advanced stage of de- 
velopment. These acoustic codes can be easily 


A Wing in Uniform Rectilinear Motion 

Consider a thin wing which is divided into 
rectangular panels and is moving parallel to the 
x^axis in uniform rectilinear motion in the x^ 2 - 
plane. Assume that the dimensions of the panel are 
2a and 2b in the chordwise and spanwise directions, 
respectively. Assume also that in the Cartesian 
frame fixed to the medium, the center of the panel 
is at (y lc , y 2c ,0). We want to the velocity 
induced at the observer position (x 1 , x 2 , x 3 ) oy a 
unit pressure distribution on the panel. Both 
integrals of the kernel function in Eg. (13) can be 
evaluated analytically since r* and t"* are explic- 
itly known. Note that 3/3n=3/3x 3 here. The normal 
derivatives of both integrals can be obtained eas- 
ily. In this example, no special treatment of the 
singularity of the second integral is needed since 
the integral is found in closed form. The alge- 
braic manipulations of this example have been car- 
ried out by Milliken [9] . 


We present the kernel function here. We first 
define the following symbols: 


e u =y lc‘V a 

V y io~ x r a 

V y 2c-V b 


V y 2c - X 2- b 


( 27-a ) 
(27-b) 
( 27-c ) 
( 27-d) 


3 


M -M 


( 27-e ) 


where M=V/c and V is the forward speed of the wing. 
The i-th panel is specified by (y 1c / y 2c ,0) and the 
control point j is specified by (x^ x 2 , 0). Ihen 
the kernel function is found to be 


K . . • — ! — ( 4ab 
13 4np 0 v n u n t 




1 r Z .2 _2 2 /72 n 2 2 . 

- — tA, +0 - /£ 0 +B n J 


( 28 ) 


- Slog [- 


s v + 


2 2 . 

V 


(Bn, + /d+6 2 n 2 )(Bn 


t + 


]} 
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where log stands for natural logarithm. 

It can be shown that Eq. (28) satisfies the 
Prandtl-Glauert rule [9] . Figure 2, shows the 
distribution of the induced velocity in the vicin- 
ity of a panel. It is very interesting to note 
that the induced velocity downstream of the panel 
is similar to that of a horseshoe vortex, as shown 
in figure 3. The pressure on a retangular wing of 
aspect ratio 10 moving at forward Mach number of 
0.2 is shown in figure 4. The Galerkin method was 
used to find this distribution. The mid-span 
chordwise pressure distribution is compared with 
the two dimensional incompressible analytic result 
in this figure. The agreement with the calculated 
pressure distribution is excellent. 

A Propeller in Uniform Axial Motion 

The kernel function for a propeller in uniform 
axial motion can be found from Eq. (23) by an ana- 
lytic-numerical method. We assume that the ob- 
server position (control point) x. never lies on 
the edge of the E-surface of a parfel. This means 
that the control point is never on the edge of a 
panel, which is a common practice in panel methods. 
Therefore, the line integrals in Eq. (23) are not 
singular. We thus concentrate on the question of 
singularities of the surface integrals only. We 
mention here that the E-surface of a panel can be 
constructed relatively easily by a numerical method 
even at transonic and supersonic speeds. This is 
done routinely in acoustic codes using time domain 
analysis [18]. 

To study the behavior of singularities of the 
surface integrals in Eq. (23) we must first set up 
a coordinate system and write the surface integrals 
more explicity. We assume that the propeller blade 
surface lies on the helicoidal surface generated by 
the blade pitch change axis (PCA). The propeller 
is moving forward at uniform speed V along the 
propeller axis and rotating at constant angular 
velocity w . The E-surface of any panel on the 
blade surface lies on this helicoidal surface. 
Essentially, one can forget about the E-surface and 
ask the following question: What is the nature of 

singularities of the surface integrals of Eq. (23) 
when we integrate over the helicoidal surface, and 
how do we evaluate these integrals? We will answer 
this question next. 

Define an n-frame such that the n 3 ~axis is 
along the propeller axis and the n^-axis is along 
the PCA at time t=0. The origin of this frame is, 
therefore, at the intersection of the PCA and the 
propeller axis^ The is defined in such a 

way that the n-frame is right handed. We first 
find the equation for the helicoidal surface. Let 
2 2 Vo 

s -n 3 and n 53 i n i + n 2 ) ' distance along 

the PCA. Since the propeller travels a distance s 
in s/V seconds, the blade rotates <us/V radians in 
this period. The helicoidal surface is, therefore, 
described by 

n=( n cos a s, nsina s, -s ) (29) 

where a=u>/V. This surface is defined by para- 
meters (n, s) • Portions of this surface are shown 
in Fig. 5. 


We next define some other geometric quantities 
on the helicoidal surface. The following vector is 
normal to the surface 


in x in 
3n 3s 


(-sina s , cos a s, an) 


its iengtn is 


.3n 3ni / 7 2 2 _ i 
|- x_| = A+a n * T( 


0(n) 


(30) 


(31 ) 


The unit normal n to the helicoidal surface is, 
therefore, 

n=0( n) (-sina s,cosa s, an) (32) 


We also have 


d E = 




dnds 

0(n) 


(33) 


Let V h be the local helical speed of the blade. 
Then and V are related by the relation 



Therefore, the local angle a^ that the helical 
velocity makes with the ^ n 2 “Pl ane is given by 

tan -1 8( n) (35) 


Assume now that the observer (or the control 
point) x. is at the distance n from the helicoidal 
surface vith the surface variables • From 

the above results we have 


Xj = n(n 0 ,s Q ) + nn(n 0 /S) 

* (n Q cosa s Q - ng^sina s Q 
n Q sina s Q + n3 Q cosa s Q , 


(36) 


s 0 + ann o 8 o ) 


where 0 =B(n 0 ) * Assuming that the source point y 
on the E-surf ace of the i-th panel has coordinates 
(n#s) on the helicoidal surface, we have 


r - x^ y = Xj- ) 

= (n n cosa s - ng sina s - ncosa s, 

( 37-a) 

n Q sina s Q + ng Q cosa s Q - nsina s, 
s - s 0+ nan 0 B 0 ) 

2 2 2 2 2 
r = n 0 + n + n + u - 2nn n cosau 

( 37-b) 


+ 2nn0 Q sinau - 2ag Q n 0 nu 


u = s 0 - s 


(37-c) 


n* r g(n), 

cos 9 = — = — - — (n 0 sma u - anu 
2 

+ ng (cosa u + a nn^) 1 
o 0 


( 37-d ) 


A ‘ 


The two surface integrals of concern here in 
the kernel function of Eg. (23) are 
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(38 - a) 

where 3/3n stands for lim^S/Sn • We remind 

readers that integrations n are over portions of the 
helicoidal surface* The upper and lower limits of 
these integrals in the variables (n#s) are indepen- 
dent of n since the influence of this latter vari- 
able is included in the line integrals of E q* (23)* 
In the following discussions, we use Eqs. (33), 

( 37-b ) and (37-d) for dl , r and cos 0 , 
respectively* 

Because we have cos 0 in the numerator of I 1 , 
we can tolerate r 2 in the denominator of the second 
term on the right of the following equation: 

3 .COS0. _ 1 3cos0 cos 3 3r , 

3n r ~ r 3n 3n 

In fact the second term yields an improper (conver- 
gent) integral* A well-known result of potential 
theory tells us that, if, x. is in the region of 
integration, then we must isolate a small region 
around x. and integrate analytically the second 

term over this small region. The result is 
2 tt [ 3r/3n ]+ • Integration over other regions then 

must he carried out separately. It is interesting 
to note that the integration with respect to n can 
be performed analytically, and only the integration 
over s must be evaluated numerically. Our work 
with I 1 is thus finished* 

For I 2 we must proceed more cautiously* We 
note that r+0 as ( r\, s ) -►( ) . We first inte- 

grate the integral with respect to n analyti- 
cally. Then we take 3/3n inside the inteqral and 
study the convergence of the resulting integrals 
We find that integrals are convergent except near 
the control point x_. . For this small region, we 

again keep 3/3n outside the integrals which involve 
now only integration with respect to the variable 
s. Using Taylor series expansion of the integrands 
in u=s-Sq, the integrals are evaluated analytically 
and then 3/3n is calculated* The result is no 
lonqer singular as n-*0. The algebraic manipula- 
tions are very lengthly and tedious but straight- 
forward. Here we have given the important equa- 
tions and the crucial steps for finding the kernel 
function for a propeller. 

Now we summarize the procedures for obtaining 
the kernel function K^<* The first two integrals 
in Eq. (23) are evaluated only once for each at 

the time t. We take this time as t^O since Ap is 
independent of time. The steps for the first two 
integrals are as follows: 

i) Construct the E-surface of the i-th panel 
numerically for t=0. 

ii) Perform the line integral over the edge 
of the E-surface numerically . 

Hi) With the surface integral written as in 
Eq, (38-a), check to see if ^ is on the E-surface. 

If so, isolate ^ by a small region and integrate 

analytically. This gives a single finite term. 
For the remainder of the E-surface, integrate the 
surface integral analytically with respect to n and 
numerically with respect to s. If x. is not on the 


E-surface, repeat the above analytic (wrt n) - 
numerical (wrt s) integration for the surface inte- 
gral. 

For the last two integrals of Eq. (23), we 
have a time integration over t" from -» to 0. We 
discretize the time integration. For each t'^<0 t 
we perform the following steps: 

i) and (ii) as above for the first two 

integrals . 

iii) write the surface integral as in Eq. 

(38-b). Integrate analytically with respect to 

n. If x_. is on the E-surface, isolate x_. by a 

small region and integrate analytically with re- 
spect to s. Take 3/3n analytically. For the 
remainder of the E-surface, take 3/3n inside and 
integrate numerically with respect to s. If x_. is 

not on the E-surface, take 3/3n inside the integral 
and integrate numerically with respect to s. 

iv) Repeat the above three steps for the next 
t and finally integrate numerically both the 

line and surface integrals with respect to t^. We 
mention here that when t'<<0 , we can simplify the 
integrands considerably and thus evaluate the line 
and surface integrals more efficiently. 

Concluding Remarks 

In this paper, we have developed the theoret- 
ical foundations necessary to solve for the aerody- 
namic loading on thin lifting bodies in arbitrary 
motion. We treat the problem from an acoustic 
viewpoint in a frame of reference fixed in the 
undisturbed medium, and we exploit the fact that 
the pressure satisfies the wave equation in this 
reference frame. Generalized function theory is 
applied to obtain an inhomogeneous wave equation in 
which the boundary conditions on the moving body 
are explicitly included as source terms. This 
equation is then solved formally using the free 
space Green’s function for the wave operator. We 
use the formal solution for the observer located on 
the body surface to derive two linear singular 
integral equations which are satisfied by the sur- 
face pressure. One of these is completely general, 
but we find that the second (simpler) one is a 
proper integral equation only in certain special 
cases. Among these are the steady problem for 
loading due to camber and the unsteady periodic 
loading problem. 

Two examples of the use of the theory as the 
basis of a numerical panel method are discussed. 
The first is a wing in steady rectilinear motion. 
In this case, detailed numerical results obtained 
elsewhere (93 are found to compare well with the 
predictions of alternate analyses. The second 
example involves a propeller in uniform axial 
motion. Here we concentrate on an analytical dis- 
cussion of the behavior of the aerodynamic kernel. 
We find that it is possible to eliminate all singu- 
lar integrals which occur, so that the current 
method promises to lead to an extremely efficient 
numerical scheme. Actual numerical implementation 
of the theory to the propeller problem is currently 
in progress. 

Ihe theory presented here has several advan- 
tages in comparison to more conventional discus- 
sions of aerodynamics. First, of course, is the 
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fact that it is applicable to arbitrary unsteady 
motions of the body so long as they are within the 
realm of linearized theory. Because it is carried 
out completely in the . time domain, there is no 
necessity for a Fourier Transform approach with its 
attendant complications of infinite series of spe- 
cial functions. Even in the case of steady recti- 
linear motion, however, the current approach yields 
results in a manner which many will find easier to 
understand than the usual development of aerody- 
namic theory. In particular, the acoustic approach 
involves the direct solution of a well-defined 
boundary value problem with no a priori modeling of 
the flow field by distributions of singularities. 
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to edge curve 


Figure 1. Geometry of edge curve of the panel and 
local plane normal to it. 
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